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AMfiUCT 

The  two-dimensional  flow  vinder  a  free  surface  with  mm-sero 
cavitation  number  Is  considered  using  a  "point  drag"  cavity  model. 
The  drag  coefficient  Is  obtained  In  closed  form  as  a  function  of 
the  cavitation  number  and  the  curvature  at  the  cavity  end  (or  nose)> 
but  not  explicitly  of  the  depth.  The  thickness -length-ratio  of 
the  cavity  Is  calculated  numerically  and  compared  with  the  case  of 
an  elliptic  cylinder  under  a  free  surface. 

The  three  dimensional  effect  on  the  cavity  Is  also  con¬ 
sidered  by  the  approximation  of  the  cavity  to  an  ellipsoid. 

These  results  should  be  of  use  In  the  estimation  of  the 
shape  of  cavities  behind  ventilated  hydrofoil  wings  and  In  the 
subsequent  calculation  of  downwash. 
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ON  FINITE  LENGTH  CAVITIES  BENEATH  A  FREE  SURFACE 
1 .  INTRODUCTION 


Fully  cavitated  hydrofoils  under  a  free  surface  at  infinite 
Proude  number  and  non-zero  cavitation  number  in  two  dimensional  flow 
are  considered.  Using  a  method  previously  used  by  Pabula  (i960)  and 
Wu  (1957)  we  compute  the  point  drag  acting  on  the  end  of  the  lin¬ 
earized  hydrofoil-cavity  system.  Thus  we  obtain  the  relations  be¬ 
tween  drag,  length,  and  thickness  of  the  cavity  for  a  given  depth 
of  foil  and  cavitation  number. 

Rapid  progress  in  the  development  of  high  speed  hydrofoil 
boats  has  led  to  a  need  for  an  understanding  of  cavity  flows  beneath 
a  free  surface.  Although  the  case  of  zero  cavitation  number  has 
been  developed  a  great  deal  in  recent  years  [e.g.,  Jakob  Auslaender 
(1961),  V.  E,  Johnson  (1958)],  the  mathematical  difficulty  in  deal¬ 
ing  with  a  doubly  connected  domain  has  retarded  the  investigation 
of  the  case  of  finite  length  cavities  (non-zero  cavitation  number) 
under  a  free  surface. 

Wu  (1957)  pointed  out  that  the  calculation  of  the  drag  from 
the  flow  field  near  the  trailing  cavity  end  singularity  in  the 
linear  theory  has  the  same  physical  basis  as  the  calculation  of 
drag  on  the  image  body  in  the  Riabouchinsky  model  or  the  determina¬ 
tion  of  the  Jet  momentum  in  the  reentrant  Jet  model.  Pabula  (July, 
i960)  used  a  two  dimensional,  linearized  model  of  a  cavity  in  an 
infinite  fluid  symmetric  with  respect  to  the  vertical  axis.  This 
model,  which  may  be  considered  as  a  linearized  Riabouchinsky  model, 
is  called  "a  point  drag  cavity  model",  and  has  the  same  singu¬ 
larities  at  the  leading  and  the  rear  ed^e  respectively.  Pabula 
(July,  i960)  calculated  the  point  drag,  and  cavity  length- thickness - 
ratio,  and  obtained  surprisingly  good  agreement  with  the  nonlinear 
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Rlabouchlnsky  model  [Plesset,  M.  S. ,  and  Shaffer,  P.  A.,  Jr.  (19^8), 
and  Byrne  Perry  (1952)].  He  also  considered  the  effect  of  the  pre¬ 
sence  of  two  symmetrically  located  solid  boxindarles  parallel  to  the 
flow  direction  with  this  model  (October  I960). 

This  report  is  an  extension  of  the  method  used  by  Pabula 
to  the  calculation  of  the  point  drag  and  the  cavity  size  in  the 
presence  of  a  free  surface. 


2.  FORMULATION  OP  PROBLEM 

Figure  la  shows  the  cavity  and  the  coordinate  system  used  In 
this  report.  Since  in  the  linearized  theory  perturbations  are  taken 
as  very  small,  the  cavity  shown  in  Figure  la  may  be  assumed  to  re¬ 
present  the  boxmdary  of  the  thin  slit  shown  in  Figure  lb.  The  line 
y  =  d  represents  the  free  s\irface  and  the  slit  at  y  =  0,  and 
0  <  X  ^  the  cavity  in  the  z  plane,  u  is  the  x  component,  and  v 
is  the  y  component  of  the  perturbation  velocity  nondimenslonallzed 
with  respect  to  the  free  stream  velocity  U^.  The  boundary  condi¬ 
tions  for  the  analytic  complex  velocity  w  are  as  follows.  Froude 
number  is  assumed  to  be  very  large  so  that  it  can  be  approximated 
to  «.  Accordingly,  u  is  equal  to  zero  on  the  free  stirface  y  =  d. 

On  the  cavity  u  =  ~  ,  where  o  is  the  cavitation  number  [see  Tulin 
0^953)].  Since  the  flow  is  considered  to  be  exactly  symmetric  with 

I 

respect  to  the  x  =  ^  axis,  u  is  symmetric  and  v  is  ant l-syimne trie 
with  respect  to  X  =  x  -  ^  .  Hence 


u(X)  =  u(-X), 
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and 

v(X)  =  -  v(-X). 

It  then  follows  that 

v(0)  =  0 

on  the  axis  X  =  0,  or  x  =  ^  .  In  view  of  these  considerations  It 

Is  necessary  to  consider  only  one  half  of  the  physical  plane 
(Figure  Ic).  The  boundary  conditions  for  the  complex  velocity 
w(z)  =  u  -  iv  In  this  plane  may  be  summarized  as  follows: 


On  the 

free  surface 

E  D 

00 

u  =  0 

On  the 

center  line 

DC>  AE 

00 

V  =  0 

On  the 

cavity 

BC>  BA 

u  =  cr/2 

At  |z| 

ss  00 

w(z)  =  0 

On  z  = 

0,  w(z)  behaves 

like  (R.  T. 

Jones,  i960) 

On  pt. 

D 

w(z)  =  0 

There  Is  no  singularity  except  at  z  =  0  in  the  half  plane. 


3.  TKANSFORmTIOK  OF  THE  PHYSICAL  PLANK 

We  now  transform  the  forward  half  of  the  z  plane  under  the 
free  sxirface  y  =  d  Into  the  upper  half  of  the  C  plane  (Figure  Id) 
by  means  of  the  Schwarz-Chrlstofell  transformation 
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where  constants  L,  a>  and  b  will  be  determined  later.  The  boundary 
conditions  given  by  Equation  [1]  become  In  the  upper  half  of  the 
+  iTj)  plane  as  follows; 


On  -1  <  4  <  1 


u  =  {i/2 


-00  <  4  <  -1 
1  <  4  <  a 


V  =  0 


a  <  4  <  00 


u  =  0 


4  «  a 
As 

At  C  =  b 


w(C)  =  0 
w(0  -►  0 


w(4)  behaves  like 


1 

C-b 


I 


and  there  Is  no  other  singularity  except  at  4  =  b  on  the 
upper  half  of  the  4  plane. 


Notice  u  and  v  are  alternately  prescribed  on  the  real  axis. 


4.  SOLUTION 

It  Is  now  necessary  to  find  an  analytic  function  w(4)  on 
the  upper  half  of  the  4  plane  satisfying  boundary  conditions  of 
Equations  [3]  on  the  real  axis  of  4  plane.  Cheng  and  Rott  (1954) 
have  developed  a  general  method  for  solving  such  boundary  value 
problems  when  u  and  v  are  given  alternately  on  the  real  axis.  The 
first  step  Is  to  find  the  homogeneous  solution  8^(4)*  solu¬ 

tion  satisfies  the  boundary  conditions  of  Equation  [3l  with  the 
following  exceptions:  at  -1  <  4  <  1  the  condition  u  »  9/2  Is 
changed  to  u  0,  and  the  boxindary  conditions  at  ••  may  be  violated. 
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H  (O  .  lAfe*  -  D*  (C  -  a)*  , 

"o^^'  c  -  b 


u„  -  1  V- 


[4] 


where  A  Is  constant.  According  to  Cheng  and  Rott  the  uniqueness  of 
the  solution  satisfying  boundary  conditions  [3]  is  assured.  In  this 
case,  by  H^(^)  7^  0  as  but 

Ho(0  ,  , 

■  V'  >'  -►  0  as  K  00  . 


We  notice  here  from  [4] 


Vjj(?,0) 


S  -  b 


We  now  find  a  new  function 


on  -  1  <  4  <  1 

w( 


H. 


fe- 


which  Is  also  analytic  In  the 


upper  half  of  the  C  plane  except  at  certain  points  on  the  real  axis. 
The  usefulness  of  this  function  Is  that  when  combined  with  boundary 
conditions  [3I  and  Equation  [4]  Its  Imaginary  part  Is  determined  on 
the  whole  real  axis. 


w(C)  u  -  lY 

W'"h- 


[51 


According  to  the  definition  of  H^(Oj  ^  portions  where 

u(|,0)  Is  prescribed.  Thus,  equating  the  Imaginary  part  of  Equation 
[5] A  we  have 
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Im 


Pi 


-  (4-b)  I 

A(«»-l)^  U-a)^ 

on  1?|<  1 


=  0 


on  ?  >  a 


Similarly  on  the  portion  where  v  (4^0)  Is  prescribed, 
Vjj(4.0)  =  0,  v(|,0)  =  0. 

We  then  have 


Now  using  the  formula  for  the  analytic  function  on  the  half 
plane  with  the  given  Imaginary  value  on  the  real  axis  [see  e.g. 
Philip  M.  Morse,  Herman  Peshbach,  Part  I  (1953)1 


Im 


w(fe,0) 


I  -  ^ 


d4 


we  have.  Inserting  the  values  [6]  and  [7]  in  Equation  [8] 
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Now  substituting  the  expression  of  H^(0  in  Equation  [4]  into 
above  equation  and  performing  the  integration  we  obtain 


W(C)  =  U(c)  -  1V(C)  -  -  --■ii-  Si. 


.1  (e-C)(«'-l)*  («-a)* 


, .  iaiCriiltili  ' _ 1 _  ,  e  -  b  i  at 


lo(C^l)i(t-a)i 

-  2  Kl 

^  (C-b)  1 

-  K 

L  va+l 

Va+l  1 

+  (C-b) 


/  at 


[10] 


where  K  is  the  first  kind  of  complete  elliptic  integral.  Near 
C  *  b  the  last  Integral  behaves  like  (C-1o)  where  the  third  kind 
of  elliptic  integral  II  has  at  most  logarithmic  singularities. 

When  we  use  the  PlemelJ  formula  [I.  I.  Privalov  (1919)^  or 
Muskhelichvili  (1946)] 

C -"t  J 

-1  -i 


where  t  is  on  the  real  axis  (>1  <  t  <  l)  and  P  /  is  the  principal  part 
of  the  integral^  we  can  easily  show  that  Equation  [10]  Is  our  solu¬ 
tion  satisfying  all  the  given  conditions. 
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5.  PflAQ 

\ 

Using  the  formula  for  the  point  drag  by  Wu  (1957)  or  by  the 


Blaslus  theorem  vte  obtain 


^  w* 
z=0 


(z)  dz 


-  Im 


w*(C)  dC 


=  Im 


5*b 

r  l^(C»-l)*(t-a)^K  (-)^) 


?=b 


Tf^  (C-b)  (a+1) 


dC 


where  +  — - —  Indicates  the  terms  which  do  not  contribute  to  the 
residue.  Kow  If  we  take  the  Integral  along  the  small  semi-circle  at 
C  =  b  using  the  value  of  the  residue,  we  obtain 


D  r  (a+1) 


[111 


constants  a,  b,  and  L  are  determined  from  the  correspondence  between 
the  z  and  C  planes.  When  the  depth  of  the  foil  becomes  Infinite, 

we  can  show  that  as  a  -» •»,  L  behaves  like  '|’'\/a~  ;  as  b  -«•  0 

if  -*  f  *  and 

Oj-fi.*.  [121 


which  agrees  with  Tabula's  result  for  the  case  of  Infinite  submergenoe 
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Thls  la  also  shown  In  Figure  2  where  the  value  of  at  d/I  -1.5 

Is  already  close  to  ir/8. 


6.  MAXIMUM  CATITY  THiCKMBSa 


The  maxlmun  thickness  of  the  cavity  Is  obtained  by  Integrat¬ 
ing  the  y  component  velocity  on  the  cavity  boundaz^  (O  <  *  <  ^  * 
y  =  0) 


or 


t  1  f  „ 


11m 
€-►  0 


l-€ 


1-b 

i(a-l)* 


1-hl 

r  h-b 

L(a-h)4 

.  1+b 

2^/r  log 

(a+l)*  ) 

hH-b 

(a4h)^ 


dh 


~  2  (log  2)V^ 
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If  we  put  h  =  cos  0,  (h®-l)^  in  the  denominator  will  be  cancelled 
and  only  a  logarithmic  singularity^  which  is  Integrable,  remains. 
This  Integral  has  been  evaluated  numerically  on  an  IBM  1620  for 
sufficiently  small  €.  The  curve  of  t/(/o)  versus  d/i  is  shown  in 
Figure  3>  where  d  is  given  by  Equation  [l6a].  The  error  is  less 
than 


2^JT  log  €  -  k^fe  -  2  (log 


1-b  ^ 


l+b 


\(a-l)*  (a+1) 


I 


l-€-b 


irV^-e  \(a-l+e)^  (a+l-e) 


l-€+b 


As  far  as  the  approximate  value  of  the  maximum  thickness  is 
concerned,  it  may  be  obtained  easily.  We  know  from  linearized  aero¬ 
foil  theory  that  a  thin  elliptic  cylinder  in  a  uniform  flow  has  a 
uniform  pressure  on  its  surface  [R.  T.  Jones  (i960)].  Even  in  the 
presence  of  a  free  surface,  when  the  depth  la  large  the  effect  of 
the  free  surface  image  above  is  very  small.  When  the  depth  is  nearly 
zeroj  the  velocity  near  the  body  is  nearly  uniform  since  the  image 
soiurce  distribution  is  of  same  magnitude  and  opposite  sign  as  that 
under  the  free  surface.  Therefore  we  can  say  that  the  pressure  on 
the  surface  of  the  submerged  elliptic  cylinder  near  the  free  sur¬ 
face  is  still  approximately  uniform. 

I^om  the  table  given  by  Jones  (i960,  p.l8)  for  a  thin  ellip¬ 
tic  cylinder  whose  major  axis  is  at  -1  ^  x  ^  1 


u  - 


z 


where  z  =  x  +  yi 


a 


«0 


or 


2u.2| 
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where  a  denotes  the  cavitation  number  In  the  Infinite  medium,  for 

to 

the  effect  of  the  free  surface  we  have  only  to  consider  a  biplane 
Image  above  (at  y  =  d).  (O  <  x  ^  i  the  distribution  of  sources  Is 
the  same  as  that  of  the  cavity  In  magnitude  but  of  opposite  sign. ) 
Approximately  the  value  of  u  at  x  »  o. 


t  4  d/i 


f4dl 

Ml  +  1 

i  i  J 

r 

where  a. Is  the  cavitation  number  for  the  depth  d.  If  we  approximate 
a 

the  cavity  \xnder  a  free  sxirface  by  an  elliptic  cylinder 


[141 


From  this,  a  curve  of  versus  4-  Is  shown  In  Figure  3  with  the 

* 

d 


case  of  the  point  drag  model.  This  Indicates  that  the  values  of 
j  almost  coincide  In  both  cases  except  between  the  values  of  d/io.i^ 


and;. 5  where  the  point  drag  model  differs  by  at  most  5.15  from  the 
elliptic  model. 


7.  CORRBSPONDENCE  BBWaOf  z  AND  C  PLAIBB 

If  we  Integrate  Equation  [2]  from  ^  ~  ^  any  point 
0  <  X  <  j  for  y  «  0  we  obtain  from  an  Integral  table  [Wolfgang 
Qrflbner  and  Nikolaus  Hofreiter  (193^)1 
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X  -  ■!  (a-b)  P  (♦,k)  -  (a+1)  E  (♦,k);>  +  | 


r 


a+l 


[15l 


where 


^  ^  ***  <  ^  >  ^  *  i  ^  ♦  -  1  and  P  and  E 


are  incomplete  elliptic  integrala. 

We  note  the  correspondence  of  points  between  the  z  and  ^  plane 


z  =  -T- 


z  = 


I  +  id 


z  =  0 


^  =  -  1  and  1 


4  =  a 


C  =  b 


Inserting  these  conditions  in  Equation  [15]  we  obtain 

-  (a+l)  E 


0  =  -~=r  !  (a-b)  K 


a+l 


■  CH 

I  v'  a+l  / 


a+l' 


[l6a] 


d  = 


2L 


a+l 


(1+b)  Kh/g^l  +  (a+l)  Er 


i  =  - 


4l 


r 


a+l 


(a-b)  P  ♦, 


a+l  1 1 


all  -  (“«)  ®  K 


[l6b] 


[l6o] 


where 


♦  =  sin 


-1  -  /  b+1 

V  2 


K  and  E  are  complete  elliptic  Integrals. 
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We  obtain  the  relation  between  a  and  b  from  Equation  [l6a] 


Combining  Equations  [l6b]  and  [l6c] 


d 

‘Viil 

+  (a+l)  E 

-a/sEII 

V  a+l  . 

Jt  ~ 

-2 

(a-b)  P 

-  (a+l)  E 

-  cul 

Dividing  Equation  [11]  by  Equation  [l6c]  we  obtain 


[17a] 


[17b] 


(b 

-a)  K® 

1-  /ml 

l  '  /a+l  1 

-4ir 

r 

[(a-b)  P 

^ 

i*'  'y^  I 

1  -  (a+l)  ! 

/a+l 

Cjy/a^t  versus  d/i  Is  plotted  In  Figure  2.  It  shows  that 

becomes  nearly  the  same  as  the  value  at  Infinite  depth  when  the 
depth-length  ratio  of  the  cavity  (d/i)  becomes  1. 


[18] 


8.  CURVATURE  AT  THE  END  POINTS 


The  curvature  at  x  =  0  (or  ^  =  b)  la  obtained  from 


dx  1  +  0/2 
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Ihus 


C  = 


11m 
4-  b 


.  a(t^-l)  (fe-a) 
2TL(4-b)»(l+  |) 


-1 


(t-b)  dt 


1  + 


(1+  |)r(4-b) 


(t-b)  I  dt 


■^-1 


C 


^  (I  +  §f 
L(l-b^ )^(a-b)^ 


since 


dt  ^  -2  ['2 

(t*-l)^(t-a)^  '/a+1  ' 


[191 


Combining  Sqiiatlons  [19]  and  [11]  we  obtain 


(1  I)- 


[20] 


That  is,  the  drag  coefficient  is  proportional  to  the  radius  of 
curvature  at  the  nose,  for  a  given  cavitation  number.  It  is  in¬ 
teresting  to  note  that  this  relation  is  Independent  of  the  depth. 
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9.  WBIAIgS 

Since  It  was  not  necessary  to  specify  the  geometry  of  the 
body  generating  the  point  drag  cavity,  there  are  not  enou^  condi¬ 
tions  to  obtain  the  length  of  the  cavity  or  the  drag  separately. 
However,  when  the  drag  of  the  body  can  be  estimated  by  the  case  of 
Infinite  depth  (M,  P.  Tulin,  ?953)  or  zero  cavitation  number  (V.  B. 
Johnson,  195B,  or  Jakob  Auslaender,  1961)  since  drag  Is  not  sensi¬ 
tive  to  these  parameters.  It  Is  possible  to  determine  the  cavity 
length,  maximum  thickness  and  hence  the  flow  field  about  the  cavity. 


10.  THREK  DIKSBSIOMAL  EFFECT  ON  THB  CAVITY  PLOW 

Since  a  two  dimensional  cavity  can  be  approximated  by  the 
known  solution  of  the  flow  past  an  elliptic  body.  It  may  be  worth 
while  to  see  the  three  dimensional  effect  on  the  cavity  flow  by 
the  approximation  to  the  flow  past  an  ellipsoid.  Not  only  the  exact 
solution  for  the  flow  past  the  ellipsoid  (Lamb,  19^5)  but  also  the 
linear  solution  (Weber,  1957)  has  been  obtained.  For  a  foil  of 
small  aspect  ratio,  Weber  obtained  constant  velocity  distribution 
on  the  elliptic  foil  by  the  application  of  slender  body  theory.  In 
Figure  6  she  compared  the  velocity  Incitements  at  the  center  of  the 
ellipsoids  by  four  methods:  the  exact  theory,  the  slender  body 
theory,  the  linearized  slender  body  theory,  and  the  ordinary  linear 
theory  respectively.  She  also  suggested  the  better  approximation 
by  the  application  of  a  factor  to  the  velocity  on  the  ellipsoid 
obtained  by  the  linear  theory.  We  use  tiie  right  handed  xyz- 
rectangular  coordinate  where  x  Is  In  the  direction  of  the  uniform 
stream  and  z  vertically  upward.  We  consider  an  ellipsoid  whose 
three  axes  are  In  the  directions  of  xyz  axes  respectively.  Then 
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the  factor  Weber  suggested  Is  written  as 


1 


'  1  +  (dz/^x) 


This  factor  reduces  the  approximate  velocity  on  the  elliptical  cylin¬ 
der  to  the  exact  one  In  the  two  dimensional  case.  In  the  present 
paper  this  factor  is  applied  to  the  approximate  velocity  at  the 
center  of  the  ellipsoid,  i.e. 

I 

U  =  U  <1  +  (dz/^x)  - 
max  ]  '  '  '  : 

L 

If  we  approximate  the  cavity  shape  to  this  ellipsoid  in  the 
uniform  flow,  the  drag  due  to  the  cavity  can  be  obtained  by  the  In¬ 
tegration  of  the  X  component  of  the  thrust  over  the  front  half  sur¬ 
face  of  the  ellipsoid.  Prom  Bernoulli's  equation 


where  and  P  are  the  velocity  and  pressure  at  Infinity  respec¬ 
tively,  and  p  Is  the  local  pressvire  on  the  surface  of  the  ellipsoid. 

Let  the  angle  between  the  tangent  plane  to  the  ellipsoid  and 
the  X  axis  be  denoted  by  a.  'Rien 


ds 

j 

s 

1 


s 


a 


L 


\ 


u‘ 


max 


1  +  (dz/dx)*^  J 
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where  Ar  is  a  given  area,  is  the  half  front  surface  of  the  ellip¬ 
soid,  and  s  is  the  projection  of  the  ellipsoid  on  the  vertical  een- 
2 

ter  plane . 

Since  the  ellipsoid  can  be  represented  as 
z  (x^y)  =  I  '  2x)*  -  (y/s)® 

where  t  is  the  maxionun  thickness  of  the  ellipsoid,  s  is  the  half 
span,  and  length  i  is  considez^d  to  be  1. 


1  „  _ 4  z®/t® _ 

1  +  (Bz/^x)*  t®  +  4  (l-t®)  z*/t®  -  t®  (y/s)* 


The  integration  of  this  expression  on  the  projection  of  ellipsoid 
on  the  plane  x  =  0.5  is 


(i-(y/s)P 


. .  •  ■  ■  dzdy 

l+(Bz/Bx)® 


tan 


I 


Vi^ 
— ir- 


sv 

T 


Hence  if  we  represent  the  cavitation  number 

P  -  P. 


a  = 


min 


«o 


-  1 


or 


U® 

max 

U  * 


1  +  0 
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we  obtain 


^  J 

r  (1+0 )  t  -1  1 

t®  +  0  1 

Ar 

\  ^  1 

1  -  t®J 

taking  Ar  as  the  maximum  cross  sectional  area  Ar  =  In  the 

dimensional  quantity. 

Bils  shows  that  the  drag  coefficient  does  not  Include  the 
semi-span  s  explicitly  as  long  as  we  nondlmenalonallze  It  by  the 
maximum  cross  sectional  area  and  express  In  terms  of  t  and  a, 

but  a  Is  a  function  of  s  and  t.  If  we  consider  a  as  an  Independent 
variable,  t  can  be  represented  as  a  function  of  o  and  s. 

To  obtain  the  free  surface  effect  on  the  three  dimensional 
cavity,  we  may  approximate  the  three  dimensional  cavity  for  mathe¬ 
matical  simplicity  as  like  a  rectangular  wing  of  elliptic  section. 
As  far  as  the  linear  theory  Is  concerned,  the  difference  In  the 
maximum  velocities  on  the  ellipsoid  and  on  the  rectangular  wing  of 
elliptic  section  Is  very  small  (see  Weber,  1957)*  Due  to  Weber 
(1957),  the  Increment  of  velocity  due  to  the  rectangular  elliptic 

wing  In  the  uniform  flow  can  be  approximated  by  the  velocity  pro¬ 

duced  by  the  source  distribution 

q(x' )  -  au^t  -r-i.  ^  In  -s  <  y  <  s 

Jl  -  (l-2x')»[  0  <  X  <  1 

I  J  z  «  0 


l.e.  the  X  component  of  the  Increment  of  the  velocity  at  (x,y,z). 
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Vjj(x,y,z) 


0  -s 


X  -  X' 


(x-x'  )®+(y-y'  )*+z* 


dx ' dy ' 


Hence 


v^(0.5,0,0)  ^ 


=  —  ts 

T 


'q  a-(l-2x' )  >  <(1-2x' )®+4s® 


4  ts 


(l+4a®)' 


K  (k*) 


where  K(k®)  Is  the  complete  elliptic  integral  of  the  first  kind  with 


k*  = 


l+4s* 


As  in  the  two  dimensional  case,  the  image  source  distribution 
for  each  x  above  the  free  surface  is  of  the  same  magnitude  but  of 
the  opposite  sign  to  that  of  the  real  source  distribution  of  the 
wing  iinder  the  free  surface.  Hence  the  increment  of  the  velocity 
due  to  the  image  also  changes  sign,  l.e. 
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v^(0.5,0,2d) 

r  (l-2x')®  dx' 

^0  ) 
( 

j  -|l6d*+(l-2x' )®|><jl-(l-2x' ^ 

r< 

^(l-2x'  )“+l6d*+4s“ 

4ts 

r 


dx' 


<jl-(l-2x')*  \  <J(l-2x' )®+l6d*+48® 
1 


+  iS2-l6cl» 


dx' 


"  /  r  If  ^ 

0  <L6d®+(l-2x'  )*>  <l-(l-2x'  )®  <;(l-2x'  )®+l6d®+4s® 


I 


J  L 


putting  1  -  2x'  =  cos  0 


t/2 

4tak'  /  d0 


/  (l+k*®sln®0) 

0 


64t3d®?vk' 


w/z 


d0 


'b 


(l+Xsln*0)(l+k'®sln*e) 


I 


K  (k'®)  -  tsd®7^k'  n  (x,k'®) 


with 


k'  = 


_ X _ 

(l+l6d*+4s® 


and  Tv  =  - 


1  +  l6d® 

Hence  the  total  Increment  of  the  velocity  at  x  =  0.5#  y  =  0,  z  =  0  Is 
v^(0.5,0,0) 


U. 


k  K  (k*)  -  k'  K  (k'*)  -  l6d*Ak'  II  (X,k'*) 


If  vie  use  the  linear  relation 


2v^(0.5#0,0) 


U. 
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we  obtain 


t 

a 


8a/kK(k®)-k'K(k'*)-l6d*Xk'n(X,k'*) 


This  Is  computed  by  using  the  tables  of  Herbert  Bristol  Dwight  (1958) 
and  Heuman  (1941).  Figure  4  shows  the  maximum  cavity  thickness  ver¬ 
sus  depth  for  each  span.  As  In  the  two  dimensional  case,  the  maximum 
cavity  thickness,  t/(ai)  at  d/i  «  i  is  almost  the  same  as  that  at 
d/i  =  «  for  each  span.  Figure  5  shows  the  maximum  thickness  versus 
the  span  for  each  depth.  The  maximum  cavity  thickness  t/(oi)  at 
s/i  =  1  Is  almost  the  same  as  that  at  s/i  »  or  as  In  the  two 
dimensional  case  for  each  depth. 
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